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Abstract 

BL Lacertae objects are an extreme subclass of AGNs showing rapid and large-amplitude variability, high 
and variable polarization, and core-dominated radio emissions. If a strong beaming effect is the cause of the 
extreme observation properties, one would expect that these properties would be correlated with each other. 
Based on the relativistic beaming model, relationships between the polarization and the magnitude variation 
in brightness, as well as the core-dominance parameter are derived and used statistically to compare with 
the observational data of a BL Lacertae object sample. The statistical results are consistent with these 
correlations, which suggests that the polarization, the variation, and the core-dominance parameter are 
possible indications of the beaming effect. 

Key words: BL Lacertae Objects: general - galaxies: active-galaxies:jets -galaxies:nuclei-galaxies: 
photometry- Polarization 



1. Introduction 

BL Lacertae objects are generally described as a sub- 
class of active galactic nuclei (AGNs), showing rapid and 
large-amplitude variation, variable and high polarization, 
a core-dominated non-thermal continuum with some BL 
Lacertae objects showing superluminal motion and high- 
energy gamma-ray emissions ( Angel, Stockman 1980; 
Antonucci, Ulvestad 1985; Fan et al. 1996; Cheng et 
al. 2000; Ghisellini et al. 1993; Hartman et al. 1999; 
Luna et al. 1993; Romero et al. 1995; Shakhovsky, Efi- 
mov 1999; Stickel et al. 1993; Takalo 1994; Wills et al. 
1992; Xie et al. 1994). There are two subclasses of BL 
Lacertae objects, namely, the radio-selected BL Lacertae 
objects (RBLs) and the X-ray selected BL Lacertae ob- 
jects (XBLs) from a survey or high-frequency peaked BL 
Lacertae objects (HBLs) and low-frequency peaked BL 
Lacertae objects (LBLs) from the spectral energy distri- 
bution (SED). RBLs correspond to LBLs, while XBLs 
generally correspond to HBLs. 

The observational properties of RBLs are systemati- 
cally different from those of XBLs. The latter have a flat- 
ter spectral energy distribution from the radio through 
X-ray regions, a higher starlight fraction, a higher ob- 
served peak of the emitted power from the radio through 
the X-ray spectral energy distribution and convex optical 
to X-ray continua. Furthermore, XBLs show lower po- 
larization compared with RBLs, and they both occupy 
different regions in the effective spectral index plot (see 
Morris et al. 1990; Junnuzi et al. 1994; Giommi et al. 



1995; Sambruna et al. 1996; Fan, Xie 1996; Fan et al. 
1997). 

Observations show that the polarization in Mkn 421 is 
correlated with the brightness of the source (Tosti et al. 
1998); a similar phenomenon was also observed from 3C 
345 (Smith 1996). Do the observations mean that there is 
a correlation between the polarization and the variations? 
In this paper, we discuss this correlation, and explain it 
in terms of a relativistic beaming effect. In section 2, 
a relation between the polarization and the variation is 
derived, and a comparison between the prediction and 
the observed data is presented. In section 3, we give 
some discussion and a conclusion. 

2. Model 

From the work of Padovani and Urry (1990) (see also 
Urry and Padovani, 1995), the observed flux, S'^^, of a 
relativistic jet is related to its intrinsic flux, 5™, by 5°*^ — 
5^5™, where S, the Dopplcr factor of the jet, is defined 
hy 6 = [r(l — /3 cos 0)]~^, /3 is the velocity in units of the 
speed of the light, F = (1 — /3^)~^/^ is the Lorentz factor, 
and 9 is the viewing angle. The value of p depends on the 
shape of the emitted spectrum and the detailed physics of 
the jet (Lind, Blandford 1985), p = 3 -I- a is for a moving 
sphere and p ~ 2 -I- a is for the case of a continuous 
jet, where a is the spectral index. However, because 
the emissions of an AGN can not be from the jet alone, 
it is natural for one to suppose that they are from two 
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components, namely the beamed and the unbeamed ones. 
In this model, the observed total flux, S°^ , is the sum of 
unbeamed, Sunh and beamed, S?^ parts. Following the 
work of Urry and Shafer 1984 (see also Urry, Padovani 
1995), we assume that the intrinsic flux of the jet, S'j'^, is 
some fixed fraction, /, of the unbeamed flux, S'unb, i-e., 
5*]" = fSunh-, we have 3°'' = (1 + fSP)S^nh- For the 
case that / is not a constant, its effect on our results 
is discussed in section 3. Furthermore, the flux is not 
totally polarized in the jet, and it is not unreasonable to 
assume that the jet flux consists of two parts, namely, 
the polarized and unpolarized parts, with the two parts 
being proportional to each other, i.e., S*]" = Sjp + 5'jup, 
Sjp = ?75'jup, where 77 is a coefficient which determines 
the polarization of the emission in the jet. The observed 
optical polarization can thus be expressed as 



pob __ 
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where the intrinsic polarization P'", is defined by 

l+fl+V 
and So is the optical Doppler factor (Fan et al. 1997) 



(2) 



2.1. Correlation between the Polarization and the Vari- 
ation 
Using 5'°'^ = (l+/(5P)S'u„b = S'olO"°''"°'' and relation 
(1), we can obtain 
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from which follows 
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for the observations of any two epochs, where P°^ and 
771°'' are the observed polarization and magnitude, re- 
spectively, while fc is a constant of proportionality. Us- 
ing S°^ = (53+"S''" + S'unb, the ratio ( = {^f+°' can be 
expressed as 
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if Si is greater than 5*2 • In this sense, the ratio C can 

gob 

be written in the form ( — (-^) , where the parameter A 

can be expected to be near unity, since (5'^+" 5'" is usually 
much greater than S'unb- Thus, the ratio of polarizations 
(4) yields 
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where Am is the variation of magnitude, which follows 



log Pi == 0.4(A - l)Am + log P2 



(6) 



Relation (5) or (6), in principle, can be investigated 
based on the simultaneous observations of polarizations 
and magnitudes (or flux density). Unfortunately, for 
most objects, there are no simultaneous observations for 
polarization and magnitude ( or flux densities). In this 
case, if we set P2 to be the minimum polarization, Pmin, 
we obtain 



log P(%) = 0.4(A - 1)Ato + log P^ 



(7) 



In order to avoid a possible observational bias (see sec- 
tion 3), we adopt the median polarization (PMod) as P 
and the largest amplitude variation (A?7iMax) as the vari- 
ation Am in our discussion. The half value of the sum of 
the maximum and minimum polarization is taken as the 
median polarization (Pmcci)- 

2.2. Results 

From the available literature, we have compiled the 
corresponding maximum and the minimum optical po- 
larizations to obtain the median polarizations, largest 
amplitude variation, and the core-dominance parame- 
ters for 35 BL Laccrtae objects. They are listed in ta- 
ble 1, in which Column 1 gives the name. Columns 2 
and 3 give the largest amplitude variation and the corre- 
sponding references; Columns 4 and 5 give the maximum 
and minimum optical polarization and the correspond- 
ing references; Columns 6 and 7 give the core-dominance 
parameter and the corresponding references. The core- 
dominance parameter, R, is defined as the ratio of the 
radio core fiux density to the fiux density in the extended 
lobes, and can be expressed as a function of the Doppler 
factor ((5) and /, i.e., R = fSP (Ghisellini et al. 1993). 

For PKS 1219+285 (ON 231), which has been observed 
for about 100 years, the early data observed by Wolf 
(1916) show that the object was as bright as 12 mag- 
nitude, which results in a largest amplitude variation of 
5.4 mag. However, the 3 points discussed in the paper by 
Wolf (1916) are perhaps not certain, since they deviate 
from other observations by about 2.3 mag, if we do not 
take the 3 early points into account. The largest varia- 
tion is thus only 3.13 (Fan, Lin 2000a), which is adopted 
in the present work. For 1400+162, Jiang et al. (1999) 
obtained the VLBI total (165 mJy) and core (114 mJy) 
fluxes at 5 GHz, suggesting a core-dominance parame- 
ter of i? — iQ^^ii4^ = 2.2. The relevant points are 
shown in flgure 1 . When a least-square regression fitting 
is performed to the median polarization and the largest 
amplitude variation for these 35 objects, the obtained 
result is 



log P{%) = (0.10 ± 0.01)Am + 0.72 ± 0.04 



(8) 



No. 2] 



Polarization and Variations of BL Lacs 



Table 1. Observation data of BL Lacertae Objects 



Na me 


Am-Max 


Ref 


-PMax 


Ref 


logi? 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


0048-097 


2.7 


FL99 


27.1/7.0 


AS80, CH84 


0.97 


W92 


0118-272 


1.05 


FL99 


18.7 


SF97 






0215+015 


5.0 


P83 


20.0/11.9 


AS80, B86 


0.90 


W92 


0219-1-428 


2.97 


FLOOB 


18.0/6.0 


FL99, RS85 


0.25 


W92 


0235-^164 


5.3 


FL99 


43.9/6.0 


ST93, FL96 


2.25 


W92 


0323+022 


1.3 


F86 


10.4/1.0 


J94 






0521-365 


1.4 


AS80 


11.0/2.8 


B83, B86 


0.01 


W92 


0537-441 


5.4 


FLOOA 


18.8/10.1 


IT90, IT88 


2.3 


G93 


0716+714 


5.0 


F97 


29.0/13.9 


ST93, RS85 


0.88 


T99 


0735+178 


4.6 


FL99 


36.0/3.0 


IT90, FL96 


3.4 


W92 


0754+100 


3.16 


FLOOA 


26.0/4.0 


IT90, FLOOA 


1.14 


W92 


0818-128 


3.78 


FLOOA 


36.0/8.0 


FL99, FLOOA 


0.23 


W92 


0823+033 


1.41 


FLOOA 


23.7/3.7 


IT88, KS90 


0.8 


G93 


0828+493 


2.0 


B90 


7.9/1.4 


KS90 






0829+046 


3.58 


FLOOA 


20.5/12.0 


VW99, CH84 


1.07 


W92 


0851+202 


6.0 


F98A 


37.2/1.0 


ST93, AS80 


3.5 


W92 


1101+384 


4.6 


FL99 


16.0/0.0 


T98, RS85 


1.0 


G93 


1144-397 


1.92 


FL99 


8.5/0.0 


IT88 






1147+245 


0.99 


FLOOA 


13.0/3.0 


IM82, IT88 


1.42 


W92 


1215+303 


3.1 


FL99 


14.0/8.0 


W78, W92 


0.27 


W92 


1219+285 


3.13 


FLOOA 


20.0/2.0 


E99, RS85 


3.45 


W92 


1308+326 


4.17 


FLOOA 


28.0/0.0 


ST93, MS81 


1.70 


W92 


1400+162 


2.8 


Z81 


14.0/4.0 


AS80, CH84 


0.34 


ST 


1418+546 


4.8 


FL99 


24.0/2.0 


FL99, AS80 


1.77 


W92 


1514-241 


2.5 


FLOOA 


8.0/2.0 


ST93, AS80 


2.41 


W92 


1519-273 


2.43 


FL99 


11.4/5.9 


IT90, IT88 


0.9 


G93 


1538+149 


3.7 


FL99 


32.8/5.0 


B86, M90 


0.95 


W92 


1652+398 


1.3 


FL99 


7.0/2.0 


CH84, KS90 


1.8 


G93 


1727+502 


2.1 


FL99 


6.0/2.0 


CH84 


1.01 


W92 


1749+096 


2.7 


FL99 


32.0/3.0 


B86, KS90 


2.83 


W92 


1749+701 


1.40 


B90 


20.3/3.5 


W92, RS85 


1.10 


W92 


1807+698 


2.23 


FLOOA 


12.0/0.0 


AS80, KS90 


0.6 


G93 


2155-304 


1.85 


FLOOB 


14.2/2.2 


PE96, B86 


0.66 


W92 


2200+420 


5.31 


F98B 


23.0/2.0 


ST93, AS80 


2.41 


W92 


2254+074 


3.07 


FLOOA 


21.0/8.8 


AS80, KS90 


1.77 


W92 



AS80: Angel, Stockman (1980); B83: Bailey et al. (1983); B86: Brindle et al. (1986); B90: Bozyan et al. (1990); CH84: 
Cruz-Gonzalez & Huchra (1984); E99: Efimov (1999); F86: Feigelson et al. (1986); F97: Fan et al. (1997); F98A: Fan et al. 
(1998a); F98B: Fan et al. (1998b); FL96: Fan, Lin (1996); FL99: Fan, Lin (1999); FLOOA: Fan, Lin (2000a); FLOOB: Fan, Lin 
(2000b); G93: Ghisellini et al. (1993); IM82: Impey et al. (1982); IT88: Impey, Tapia (1988); IT90: Impey, Tapia (1990); 
J94: Jannuzi et al. (1994); KS90: Kuhr, Schmidt (1990); MS81: Moore, Stockman (1981); M90: Mead et al. (1990); P83: 
Pettini et al. (1983); PE96: Pesce et al. (1996); SF97: Scarpa, Falomo (1997); ST: see the text; ST93: Stickel et al. (1993); 
T99: Tian et al. (1999); T98: Tosti et al. (1998); VW98: Visvanathan, Wills (1998); W78: Wardle (1978); W92: Wills et al. 
(1992); Z81: Zekl et al. (1981) 
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Fig. 1.. Plot of polarization ( PMcd. (%)) against optical variation 
(Am). 



with a correlation coefRcient of r = 0.641 and a chance 
probabiUty oi p — 6.32 x 10^'"'. The best-fit result shown 
in figure 1 with a solid line implies that the parame- 
ters A = 1.25 ± 0.12 and P^in = 5.2%. In addition, 
the relation of C, gives log i?2 = 0.4AAr7i + log Ri, since 
R = fS^^"", which can be tested by observations if we 

let i?2 = -RMax, ^1 = Rmin, and Am = AmMax- 

Therefore, log i?Max = 0.4A AmMax + ci; here cl is a 
constant associated with the minimum core-dominance, 
^min- The fitting result, A = 1.25 ± 0.12, predicts that 
log i?Max — (0.50 ± 0.05)AmMax + cl. When a linear 
regression was performed on the core-dominance param- 
eter and the largest variation listed in table 1, a best-fit 
resuh, log R = (0.31 ± 0.11)Am 4- 0.39 ± 0.40, with a 
correlation coefficient r = 0.45 and a chance probability 
p ~ 1% was obtained. The best-fit result 0.31 ±0.11 does 
not conflict with the predicted slope 0.50 ± 0.05. 

3. Discussion and Conclusion 

BL Lacertae objects are characterized by the obser- 
vational properties mentioned in the introduction. The 
beaming model was adopted to explain both the par- 
ticular observational properties and some observational 
differences between RBLs and XBLs (see Xie et al. 1991; 
Fan et al. 1993, 1997; Fan, Xie 1996; Georganopoulos, 
Marscher 1999), although the viewing angle alone can not 
explain all of the difference between the two subclass BL 
Lacertae objects (Sambruna et al. 1996). To discuss the 
intrinsic properties of BL Lacertae objects, one should 
know the boosting factor. Many authors have worked on 
this topic (Xie et al. 1991; Ghisellini et al. 1993; Lah- 



teenmaki 1999; Fan et al. 1999). Recently, using the time 
variability, we determined the boosting factor and other 
physical parameters of Blazars (Cheng, et al. 1999). 

It is believed that the particularly observed properties 
of BL Lacertae objects are associated with the beaming 
effect; if so, there should be a correlation between those 
properties. In 1996, Smith reported that the polariza- 
tion in 3C 345 is strongly correlated with the brightness. 
Very recently, Tosti et al. (1998) also found that there 
is a clear correlation between the variation of polariza- 
tion degree and the brightness during the Mkn 421 1997 
outburst. A similar situation has also been found in 3C 
66 A (Efimov, Shakhoskoy 1998a) and ON 231 (Efimov, 
Shakhoskoy 1998b). 

However, for most objects, the polarization and magni- 
tude variations are not observed simultaneously. In addi- 
tion, from an observational point of view, it is known that 
bright sources are observed frequently by photometry and 
polarimetry. In this sense, there is a tendency to detect a 
larger variation and higher polarization in those brighter 
sources. This kind of bias would result in that those ob- 
jects with a larger variation are also those with a higher 
polarization. To reduce this kind of bias, we have studied 
the median polarization of these objects. We chose the 
half value of the sum of the maximum and minimum po- 
larization as the median polarization and used this value 
to investigate the correlation between the polarization 
and the variation for a sample of 35 BL Lacertae objects. 
Nevertheless, we have also used the observed maximum 
polarization to analyze equation (7) ; the correlation coef- 
ficient is 0.65. Although this chance probability is better 
than that of using the median polarization by a factor of 
2, a possible observational bias may exist. 

It is known that XBLs are not as strongly beamed, and 
that the observed data can be taken as the intrinsic data 
to some extent (Fan, Xie 1996). In this sense, the polar- 
ization of XBLs, which is ^ 5% on average (see Jannuzi 
et al. 1994; Fan 1999), can be taken as the minimum 
polarization of BL Lacertae objects. The best-fit result 
for the polarization and the variation implies a minimum 
polarization of 5.2%, which is consistent with the obser- 
vation result of XBLs (Jannuzi et al. 1994). Our results 
indicate that the polarization is correlated with the varia- 
tion, which is also consistent with the observation results 
in both 3C 345 and Mkn 421. 

The polarization is found to be associated with the 
core-dominance parameter (see Wills et al. 1992 and 
reference therein) with a high polarization corresponding 
to a large log R. From relation (3), one can obtain a 
relation between the polarization and the core-dominance 
parameter for a certain magnitude: 



P°^ = k (5^+" 10° 
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Fig. 2.. Plot of polarization (PMax%) vs core-dominance param- 
eter (log R) 



where c{m) — (4)10"'^™ is a parameter that depends 
on the magnitude. The relation shows that the high po- 
larization is associated with a large core-diminance pa- 
rameter. 

From the adopted c{in) — 0.01, 0.1, 1, 10, several 
curves are obtained (figure 2); the curves fit the observa- 
tion data well. The differences in c{m) are from the dif- 
ferences in the magnitudes amongst the objects. From a 
catalogue by Hewitt and Burbidge (1993), the maximum 
magnitude difference among the considered objects is 5 
mag, which gives a difference of 100 in c{m); this value 
does not conflict with the difference of c{m) adopted in 
the present paper, because BL Lacertae objects are vari- 
able. The variability of the objects, themselves, should 
result in a larger than 100 difference in c{m). In addi- 
tion, although we used a fixed / in our discussion, / is 
likely to be variable, which would result in a dilution of 
the correlation (cf. figure 2). Besides, the variation in rj 
would also affect the correlation. 

We have tried to separate the data into three sub- 
groups, i.e., the lower (Am < 2.0 mag), intermediate ( 
2.0 mag < Am < 4.0 mag), and higher (Am > 4.0 
mag) amplitude variation groups, and found that the 
higher amplitude variation group shows a better linear 
relation between the polarization (P) and the variation 
(Am). We suspect that because the P for the sources in 
this group may be obtained when the sources are bright, 
they have similar magnitudes {m). However, the sample 
of such a group is small and the statistical significance is 
very low. More observations in the bright state should 
confirm the result. 

In this work, the correlation between the polarization 



and the variation was derived and explained in the rel- 
ativistic beaming frame for a BL Lacertae object sam- 
ple. Our statistical results show that the particularly 
observed properties of BL Lacertae objects, such as the 
large-amplitude variation, core-dominance, and high po- 
larization arc associated with the relativistic beaming 
model. The mutual relationships are consistent with 
the beaming model, suggesting that those observational 
properties are possible indications of a beaming effect. 
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